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Abstract 


The  time-dependent  behavior  of  a  transversely  magnetized,  twodimensional 
plasma-wall  sheath  has  been  studied  through  particle  simulations,  which 
have  shown  that  the  cross-field  sheath  develops  into  a  turbulent  boundary 
layer,  with  strong  potential  fluctuations  and  anomalous  particle  transport. 

The  driving  mechanism  is  the  Kelvin-Helmholtz  instability,  which  arises 
from  the  sheared  particle  drifts  created  near  the  wall.  The  sheath  acquires 
an  equilibrium  thickness  lz  ~  5  /»,,  and  maintains  large,  long-lived  vortices, 
with  amplitudes  64>  ~  -2 T,/e,  which  drift  parallel  to  the  wall  at  roughly 
half  the  ion  thermal  velocity.  The  sheath  also  maintains  a  large,  spatially- 
averaged  potential  drop  from  the  wall  to  the  plasma,  in  sharp  distinction 
with  the  unmagnetized  sheath,  where  the  plasma  potential  is  higher  than  at 
the  wall.  Accompanying  the  vortices  is  a  spectrum  of  shorter-wavelength 
fluctuations  which  induce  an  anomalous  cross-field  transport.  A  central 
simulation  result  is  that  for  >  2 uct,  the  transport  scales  like  Bohm 
diffusion,  a  result  for  which  we  have  a  qualitative  analytic  model. 


‘This  is  a  shorter  version  of  UCB/ERL  M88/2I,  20  March  1988  which  has  the  same  title.  It  has  been 
submitted  to  Physical  Review  Letters. 
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Main  Text 


The  following  report  on  our  particle  simulations  of  the  magnetized 

plasma-wall  she?+^  T'his  is  a  study  u '  plasma  transport  perpendicular  to 
a  magnetic  field,  in  a  plasma  bounded  by  a  conducting  wall,  with  the  aim 
of  modelling  the  vicinity  of  the  limiters  and  walls  of  magnetized  plasma 
devices.  Our  approach  has  been  to  use  our  two-dimensional,  bounded  par¬ 
ticle  simulation  code  ES2  to  investigate  the  edge  effects  which  arise  in  such 
a  configuration. 

Our  simulations  have  shown  that  the  cross-field  sheath  between  a  wall 
and  a  plasma  is  not  a  static  structure,  but  is  in  fact  a  turbulent  boundary 
layer,  with  strong  potential  fluctuations  and  anomalous  particle  transport. 
The  driving  mechanism  is  the  Kelvin-Helmholtz  instability  which  arises 
from  the  sheared  particle  drifts  created  near  the  wall.  Provided  it  is  re¬ 
plenished  by  an  internal  flux  of  particles,  the  sheath  will  maintain  itself  in 
a  dynamic  equilibrium,  in  which  the  linear  edge  instability,  the  nonlinear 
plasma  flows  and  the  outward  particle  diffusion  all  balance  each  other.  It 
is  important  to  emphasize  that  the  turbulent  behavior  of  the  sheath  is  a 
completely  spontaneous  phenomenon,  which  arises  from  the  self-consistent 
plasma-wall  interaction,  and  which  does  not  require  the  imposition  of  ex¬ 
ternal  fields.  This  self-consistency,  and  the  short  space  and  time  scales 
inherent  in  the  plasma-wall  interaction,  sharply  distinguish  our  work  from 
previous  simulations  of  the  Kelvin-Helmholtz  instability [l ,2,3]. 

In  Figs.(l),  we  show  a  series  of  snapshots  of  the  electrostatic  potential 
<t>{ x,  y,  t),  in  a  system  initially  filled  with  a  uniform  distribution  of  electrons 
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and  ions.  Our  simulation  is  two-dimensional  and  electrostatic!*!],  with  the 
particles  moving  in  the  (x,  y)  plane,  and  subject  to  a  magnetic  field  B  =  iB. 
Periodic  boundary  conditions  are  imposed  at  the  top  and  bottom  of  the 
simulation  region  (y  =  0,  Lv).  The  wall,  at  x  =  0,  is  a  perfect  conductor 
and  a  perfect  absorber  of  electrons  and  ions:  when  a  particle  hits  the  wall, 
its  charge  is  immediately  accounted  for  as  mobile  surface  charge.  The 
boundary  conditions  on  the  right-hand  side,  r  —  Lz,  simulate  a  semi-infinite 
inner  plasma:  x  =  L,  is  an  equipotential,  and  for  the  particles,  an  inversion 
symmetry  condition  is  applied [5] .  Finally,  there  is  also  a  distributed  plasma 
source,  by  which  warm  electron-ion  pairs  are  created  at  a  constant  temporal 
rate,  and  with  spatially  random  occurence.  For  Figs.(l),  we  have:  m,/me  = 
40,  T,  =  r„  Ljpi/Ud  =  1.15,  c Jpt/uet  =  0.18,  Pi[pt  =  6.32,  and  a  system  size 
Lz! Pi  =  5.06  and  Ly/p,  =  40.5.  Ne  —  N,  =  15000  numerical  particles 
of  each  species  were  initially  in  the  system,  with  the  creation  rate  (40  e-i 
pairs  per  u;”1  time  interval)  adjusted  to  achieve  similar  steady-state  values, 
Ne  ~  16500,  N,  «  14000. 

Fig. (la),  at  ur,t  =  15,  shows  the  essentially  y-uniform  sheath  which  has 
formed  after  just  2j  rotations  of  the  ions.  This  sheath  is  due  to  the  initial 
loss  of  ions  into  the  wall,  resulting  in  a  layer  of  depletion  of  positive  charge 
in  the  plasma  and  a  net  positive  charge  on  the  wall.  The  result  is  a  large 
potential  drop  from  the  wall  into  the  plasma,  eA<f>/T,  %  —1.1. 

The  sheath  in  Fig. (la)  contains  a  non-uniform  electric  field  Ex(x)  >  0, 
with  maximum  intensity  at  the  wall,  which  induces  a  downward  E  x  B  drift 
of  electrons  and  ions  with  a  maximum  velocity  very  close  to  the  ion  thermal 
speed,  vv(x  =  0)  =  —  Ez[x  =  0 )jB  «  —  vti.  Because  the  electric  field  is 
nonuniform,  the  resulting  flow  is  strongly  sheared,  and  is  vulnerable  to  the 
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Kelvin-Helmholtz  instability.  In  the  subsequent  evolution,  noise  grows  into 
larger  perturbations,  which  are  both  amplified  and  convected  by  the  E  x  B 
flow.  In  Fig. (lb),  the  growing  waves  have  amplitudes  e54>/T,  %  —0.2,  and 
wavelength  parallel  to  the  wall  Au  ss  60. 

As  the  growing  modes  saturate,  they  form  small  vortices  with  approxi¬ 
mate  spacing  A„  «  60  (Fig.(lc)).  These  vortices,  generated  in  60  <  uc,t  ~ 
200,  are  short-lived  because  they  rapidly  coalesce  with  each  other,  forming 
a  smaller  number  of  vortices  with  Av  ss  120.  There  is  also  simultaneous 
competition  from  the  mode  with  A„  s:  120,  which  has  a  comparable  linear 
growth  rate.  The  overall  effect  is  to  favor  the  longer  wavelength,  and  by 
uj„t  =  250,  Fig. (Id),  only  two  large  vortices,  each  with  e6<p/T es  —1  and 
Av  120,  have  survived. 

The  two-vortex  state  of  Fig. (id)  is  only  quasi-stable,  because  the  vor¬ 
tices  are  once  again  vulnerable  to  the  coalescence  instability,  and  Figs.(ld,e,f) 
(u*'eit=  250,  375  and  437.5)  show  the  progression  to  a  single,  large-amplitude 
vortex.  This  final  vortex  remains  very  stable  throughout  the  remainder  of 
the  simulation,  which  we  ended  at  u„t  —  1500.  It  drifts  parallel  to  the  wall 
at  a  near-constant  velocity  v0v  =  — 0.44v(>,  and  maintains  e6<p/T ,  %  -2.2. 
However,  the  background  plasma  remains  fluctuating,  and  from  time  to 
time  a  smaller  vortex  is  generated  at  some  distance  from  the  main  vortex, 
eventually  merging  with  the  latter. 

Overall,  we  have  explored  a  range  of  parameters  0.5  <  w’./w*  <  10, 
with  ui  l «  1,  and  have  found  behavior  qualitatively  very  similar  to 
the  one  outlined  above.  Simulations  of  broader  systems  than  in  Fig.(l)  also 
show  that  the  steady-state  vortex  and  sheath  acquire  a  “natural”  thickness 
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of  lz  ss  5  and  do  not  grow  indefinitely;  thus  the  vortices  are  a  feature  of 
the  sheath  and  not  the  bulk  plasma. 


To  explain  the  evolution  seen  in  Figs.(l),  we  have  investigated  the  fluid 
cross-field  equations.  Assuming  jw|  <  ojci,  we  approximate  the  electron 
motion  by  the  E  x  B  drift,  the  ion  motion  by  the  ExB  and  polarization 
driftsliSj,  and  obtain  the  coupled  nonlinear  equations: 


■£*’* 


V  •  n— -Vp 
1  dt 


-  C 


(1) 


(2) 


where  didt  =  d  I  dt  -4-  (E  x  B  /  B)  •  V  is  the  total  derivative  and  V  =  ( dz ,  dv) 
We  then  linearize  Eqs.(l,2),  in  the  presence  of  equilibrium  profiles  v0(x)  -- 
VoC'q(x)  and  n0(i)  =  Ar0n0(x),  with  V0  the  velocity  at  the  inflection  point 
and  No  the  bulk  plasma  density,  and  t}0(x)  and  rz0(x)  dimensionless  profiles. 
With  a  =  =  Ar0),  we  have  the  stability  equation: 


d2  ,  wp(g)  g  ,  (  ^  (  „f, 

dx1  1  a  4-  n0(x)  dx  1  -  kvv0[x)  y  0 


™c(x)vo(z) \ 
a  4-  n0(x)  ) 


Pi  =  0 

(3) 


Taking  t>0(x)  =  ta nh(x/a  -  1)  —  1,  and  n0(x)  =  tanh (x/d),  and  fitting  the 
parameters  V0,  No,  a  and  d  to  the  conditions  in  the  sheath  at  ucit  =  2 tt , 
we  find  Vo  —  0.51utl,  a  =  0.63p,,  d  -  1.26 p,  and  a  =  0.73.  We  then 
numerically  solve  Eq.(3)  for  u>  —  uR  4- 17  as  a  function  of  kv[ 2!.  The  results 
are  shown  in  Fig. (2),  where  they  are  compared  with  the  simulation  results. 
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The  predictions  for  the  real  part  of  the  frequency  agree  reasonably  well. 
On  the  other  hand,  while  the  growth-rates  for  the  long-wavelength  modes 
(m  =  1,2)  are  in  agreement,  fluid  and  simulation  results  diverge  at  the 
shorter  wavelengths  (rn  =  3,4,5,  6),  an  effect  probably  due  to  finite-larmor 
radius  effects,  but  which  we  have  not  further  investigated. 

Turning  to  the  nonlinear  regime,  we  have  attempted  to  model  the  final 
vortex  by  an  approximate  solution  of  Eqs.(l,2).  We  take  the  Navier-Stokes 
limit  of  Eqs.(l,2)  (n  =  constant,  or  >  u^,),  and  then  consider  solutions 
which  satisfy  in  a  co-moving  frame  V2d>  =  /(<£)  [3,6).  Solutions  with  linear 
f[o)  —  a  /3<p,  as  in  1 3 ,6) ,  lead  to  vortices  with  amplitudes  e6\<j>\/Tt  1, 

in  contradiction  with  simulation  results.  This  suggests  looking  for  more 
tightly-bound  solutions,  for  instance  with  an  exponential  /(<?).  In  particu¬ 
lar,  the  vortex  chain  of  Stuart[3,7]: 


d>(x,  y)  =  -v0x  +  ~  log  (cosh(fcv(x  -  b))  4-  A  cos(fcy(y  -  n0f)))  (4) 

satisfies,  in  its  rest-frame,  the  equation  V2o  =  kvv0  exp(-2kyd>/u0)-  With 
kv  -  2;r  jLv  =  0.155/pi,  vQ  —  — 0.44u(,,  b  =  2.5 p,  and  .4  =  0.57,  we  have 
found  that  Eq.(4)  gives  a  rough  but  sufficiently  large-amplitude  approxi¬ 
mation  to  the  observed  vortex  structure. 

An  essentia!  feature  of  the  cross-field  sheath  is  the  existence  of  a  nearly 
constant,  ambipolar  transport  of  particles  to  the  wall,  for  which  we  shall 
sketch  a  scenario:  we  believe  that  the  ion  motion,  with  a  gyro-radius  an 
appreciable  fraction  of  the  vortex  dimensions  ( pjlz  w  1/5),  is  intrinsically 
stochastic,  and  readily  leads  to  an  outward  ion  diffusion,  which  is  only 
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stopped  by  the  positive  charge  building  up  at  the  wall.  The  electrons  on 
the  other  hand  are  more  strongly  confined,  and  it  is  thus  their  rate  of 
transport,  induced  by  the  edge  turbulence,  which  determines  the  overall 
rate  of  transport  for  both  species. 


The  orbits  of  test  electrons  show  that  their  outward  transport  is  the  re¬ 
sult  of  the  E  x  B  motion  along  broad  orbits,  trapped  in  or  circulating  about 
the  vortices,  combined  with  a  shorter-scale  diffusive  motion,  which  we  as¬ 
cribe  to  the  the  high-frequency,  short-wavelength  fluctuations  ( ! w j  ~ 

|k;p,  ~  1)  which  we  see  in  the  simulation  power  spectrum,  and  which 
accompany  the  vortices.  Particularly  vulnerable  to  scattering  are  the  elec¬ 
trons  on  the  vortex  separatrix,  their  rapid  loss  leading  to  an  evacuated 
region  in  front  of  the  vortex.  In  the  steady-state,  we  can  define  a  time  and 
space-averaged  diffusion  coefficient  Dz  =  sl‘/2n{lz),  where  lz  is  the  sheath 
thickness  (I*  =  5pt),  s  is  the  rate  of  pair  creation  (per  unit  area,  per  unit 
time),  and  n(lz)  is  the  t/-averaged  density  in  the  plasma  interior,  assuming 
a  parabolic  n(x )  m  0  <  z  <  lz.  The  results  for  Dz  from  simulations  with 
increasing  creation  rates,  and  hence  increasing  inner  plasma  densities,  are 
shown  in  Fig.  (3),  where  Op,  is  evaluated  at  the  spatially-averaged  density. 
For  Op,  >  2wct ,  the  diffusion  coefficient  has  the  Bohm-like  dependence: 


Dz  =  0.04-^  ,  Op,  >  2wct  (5) 

We  have  verified  that  the  results  in  Fig. (3)  are  independent  of  plasma  dis¬ 
creteness,  by  varying  n\2D  in  the  range  4  <  n\2D  <  56,  thereby  confirming 
that  transport  is  not  collisional. 
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To  provide  a  derivation  of  Eq.(5),  we  first  assume  that  for  Wp,  2> 
the  edge  turbulence  is  density-independent.  This  is  consistent  with  the 
uipi  >  limit  of  Eqs.(l,2).  We  then  introduce  “universal”  parameters 
Ci,  c-i  and  c3,  with  the  vortex  spacing  lv  =  Cjp,,  the  vortex  size  /„  =  c2p., 
and  the  vortex  depth  4>maz  =  c3T',/e,  with  Ci  a  40,  c 2  ~  5  and  c3  a 
2.  Assuming  a  spectrum  <  |d>(fcv)j2  >=  A{kypx)  <  |#(0)js  >,  and  an 
wrave-e!ectron  correlation  time  rc(ky)  =  oj^1T(!:ypx),  with  the  functions  A 
and  T  independent  of  B,  we  then  follow;8’  and  find  an  estimate  which 
scales  as  in  Eq.(5),  Dz  -  c0  ( T,/eB)x  with  constant  c0  given  by:  c0  = 
(l/rr)(c2c3/c1):  /0°°  dK K‘A(K)T{K).  Assuming  the  integral  to  be  of  order 
1,  we  find  Co  =  0.02,  an  estimate  qualitatively  in  agreement  with  Eq.(5). 

We  have  also  studied  the  effects  of  finite  k j,  by  tilting  the  magnetic  field 
by  a  small  angle  8  —  (mf/m,)1/'J/2  from  the  perpendicular.  With  the  field 
tilted  parallel  to  the  wall,  the  vortices  are  suppressed,  and  the  edge  turbu¬ 
lence  and  cross-field  transport  are  greatly  reduced,  but  the  sheath  main¬ 
tains  a  large  potential  drop,  c£<t>/T,  a  —1.4.  With  the  field  tilted  towards 
the  wail,  the  morphology  is  changed  to  that  of  the  unmagnctized  sheath, 
with  electron  flow  dominant,  and  a  potential  rise  from  wall  to  p!asma;9  . 
These  results  emphasize  that  the  Keivin-Helmholtz  vortices  arise  for  flute¬ 
like  modes,  forming  only  in  regions  of  tangency  to  walls  and  limiters,  or 
perhaps  internally,  across  a  separatrix. 

Future  simulations  of  the  cross-field  sheath  should  explore  a  larger  range 
of  parameters  than  possible  in  the  present  work.  For  instance,  we  believe 
that  longer  systems  should  accomodate  many  steady-state  vortices.  We 
were  also  limited  to  <  iO  whi'e  we  would  like  to  explore  values  as 
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high  as  =  103,  representative  of  a  fusion  environment.  A  straight¬ 

forward  modification  of  our  simulation  code,  to  allow  for  an  E  x  B  electron 
mover,  should  make  these  extended  studies  feasible  at  a  reasonable  cost  in 
computer  time. 
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Figure  Captious 


Fig.(l):  Contour  plots  of  the  electrostatic  potential  <£(x,y,t)  taken  during 
the  particle  simulation.  The  wall  is  at  the  left-hand  boundary,  x  —  0. 


Fig.  (2):  Comparison  of  8uid  theory  for  the  Kelvin-Helmholtz  instability 
(full  lines),  with  the  simulation  results  (large  dots);  a)  growth  rates;  b)  real 
parts  of  the  frequency.  In  Fig. (2b),  u0  is  the  vortex  speed,  v0  =  0.^4vtt. 


Fig. (3):  Dependence  of  the  time  and  space-averaged  diffusion  coefficient  on 
the  steady-state  density. 
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